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Understanding the electronic charge distribution around oxygen vacancies in transition metal and 
rare earth oxides is a scientific challenge of considerable technological importance. We show how 
significant information about the charge distribution around vacancies in cerium oxide can be gained 
from a study of high resolution crystal structures of higher order oxides which exhibit ordering of 
oxygen vacancies. Specifically, we consider the implications of a bond valence sum analysis of 
CerOi2 and CenCbo- To illuminate our analysis we show alternative representations of the crystal 
structures in terms of orderly arrays of co-ordination defects and in terms of flourite-type modules. 
We found that in CerOi2, the excess charge resulting from removal of an oxygen atom delocalizes 
among all three triclinic Ce sites closest to the O vacancy. In CenCbo, the charge localizes on the 
next nearest neighbour Ce atoms. Our main result is that the charge prefers to distribute itself 
so that it is farthest away from the O vacancies. This contradicts the standard picture of charge 
localisation which assumes that each of the two excess electrons localises on one of the cerium ions 
nearest to the vacancy. This standard picture is assumed in most calculations based on density 
functional theory (DFT). Based on the known crystal structure of Pr 6 On, we also predict that the 
charge in CeeOn will be found in the second coordination shell of the O vacancy. Although this 
review focuses on bulk cerium oxides our approach to characterising electronic properties of oxygen 
vacancies and the physical insights gained should also be relevant to surface defects and to other 
rare earth and transition metal oxides. 
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Non-stoichiometric phases of a number of transition 
metal and rare earth oxides have become widely studied 
in the past decade because of their potential in technolog- 
ical applications. Important examples are derived from 
the "parent" materials Ti02, Zr02, Hf02, Ce02, by dop- 
ing with metal ions or reduction by removal of oxygen 
atoms. A key scientific question in all these materials is 
the structure and electronic properties of oxygen vacan- 
cies. Ceria (cerium oxide) is a technologically important 
material with applications in high temperature electro- 
chemical devices^, catalysig ^MH, o xygen gas sensor d^^HD 
and magnetic semiconductors^^. A fundamental prop- 
erty of this material relevant to all of these application 
is its oxygen storage capacity (OSC). The material can 
rapidly take up and release oxygen through a reversible 
chemical reaction. There has been a considerable multi- 
disciplinary research effort aimed at developing a fun- 
damental picture of the microscopic processes associated 
with this reversible chemical reaction. Key questions in- 
clude: 

1. What is the origin of the reversible uptake and re- 
lease of oxygen by ceria? 

2. What is the origin and mechanism of the anionic 
conduction? 

3. What is the nature, composition and geometry of 
the catalytically active sites on cerium oxide sur- 
faces? 



4. 



What is the nature of oxygen vacancies in the bulk 
solid and on surfaces? 
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5. When an oxygen atom is removed to create a va- 
cancy, what happens to the two electrons that are 
left behind? 

Similar questions are also important in other oxides. This 
review is primarily concerned with the last two questions 
and the study is restricted to the case of bulk ceria. In 
this introduction, we briefly outline the current under- 
standing of the problem with respect to the last two ques- 
tions. 

When Ce02 is reduced to the various defective phases, 
Ce02-x, according to Eq. [T] O vacancies are formed in 
the lattice structure. 

Ce0 2 ^ Ce0 2 - x + |o 2 (g), < x < 0.5 (1) 

The crystal structure adopted by any such defective 
phase, Ce02-x, is understood to be the one that pro- 
vides the most favourable energetics for the arrangement 
of all the O vacancies within the structure. In a widely 
accepted view of the microscopic description of O va- 
cancy formation and ordering in Ce02- X phases, the two 
electrons associated with a missing O atom when an O 
vacancy forms fully localize on two of the four equivalent 
Ce 4+ ions which form a tetrahedron around the vacancy 
site as shown in Eq. IpUM! 

4Ce Ce + o ^ 2Ce Ce + 2Ce , Ce + V Q + ^0 2 (g) (2) 

where we have used the Kroger- Vink notation^ so that 
the symbols have the following meanings: Cec e - a Ce 4+ 
ion on a Ce lattice site, Oo - O 2- ion on an O lattice site, 
ion on a Ce lattice site and Vq - neutral 
O vacancy site. This mechanism, which we will hence- 
forth refer to as the standard picture, is illustrated in 
the schematic in Fig. [I] The localization of an electron 
on a Ce 4+ ion converts it to the slightly larger Ce 3+ ion 
with one electron in the 4/ orbital. In the reverse pro- 
cess where a defective phase, CeC^-a, is oxidized, two 
4/ electrons from the two neighbouring Ce 3+ ion sites 
move onto the site where an O atom is incorporated and 
then delocalize into the O 2p valence band. Thus, the 
reversible processes of oxidation and reduction in Eq. [I] 
have been considered to involve the extremal states of 
the Ce 4/ electrons in which they are fully delocalized 
and fully localized respect ivel}B^. 

There have been, broadly speaking, three different ap- 
proaches for investigating this problem. Crystallogra- 
phers have attempted to establish general principles for 
O vacancy ordering in the reduced ceria phases. In 1974, 
Martin proposed the coordination defect model 18 . He 
suggested that when an O vacancy is created in the lat- 
tice of a rare-earth oxide, the local environment of the 
vacancy undergoes relaxations in such a way that the 
four nearest neighbour Ce and six nearest neighbour O 
atoms form a stable structural entity which is referred 
to as a coordination defect. The crystal structures of 
the higher rare-earth oxides are then determined and re- 
stricted by the topology of this coordination defect. In 
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FIG. 1. Schematic of the standard picture of charge redistri- 
bution following the formation of an oxygen vacancy in CeC^. 
The tetrahedron of Ce atoms (black circles) with an O atom 
at its centre (grey (orange in colour version) circle) is shown 
along with the charges on these atoms in the simple ionic pic- 
ture description of CeCb- The process of reduction shown by 
the arrow leads to a neutral O vacancy at the centre of the 
tetrahedron (empty circle) while two of the Ce ions have been 
reduced to the +3 oxidation state. 



1996, Kang and Eyring developed a different framework 
for describing vacancy ordering based on structural ele- 
ments which are derived from the conventional unit cell of 
Ce02 in a simple way^. As will be described in more de- 
tail later, there are thirteen of these structural elements, 
called modules, and they are the fundamental building 
blocks for all the crystal structures of the higher oxides 
in the Kang-Eyring fluorite-type module theory. As a 
result, it is possible to define what has been called a 
modular unit cell in which, instead of the usual atomic 
basis, the different modules form the basis. The value of 
this approach is the simplicity and elegance with which it 
provides insight into the superstructures observed in the 
reduced higher oxides in relation to the parent fluorite 
structure from which they are derived. 

Neither crystallographic approach sought to explicitly 
account for the charge redistribution that occurs when 
an O vacancy is formed. However, as we will show, they 
provide useful vantage points from which to examine the 
problem. In contrast, the remaining two approaches in 
the strategy sought to specifically establish the nature 
and occupation of the Ce 4/ level in the Ce oxides and 
consequently, the charge redistribution that occurs when 
an oxygen vacancy is created. The earliest literature in 
this direction emerged from spectroscopists^ with their 
results being interpreted using either cluster models^ or 
the single impurity Anderson modeP^HHl This work in- 
vestigated the occupation of the Ce 4/ level in the two 
extremal phases, Ce02 and Ce203 and it was concluded 
that Ce02 is mixed valence whereas Ce203 is a pure 4/ 1 
configuration. The result that Ce02 is mixed valence 
is inconsistent with the ionic description underlying the 
standard picture. However, the models used to interpret 
the spectroscopic results are not parameter- free and the 
conclusions drawn from them have been contested 25 . 
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Recently electronic structure calculations based of den- 
sity functional theory (DFT) have been performed. In 
general, DFT results for Ce02 assume empty Ce 4/ 
states which is the simple ionic picture of bonding in 
this oxide, e.g., RefPP. To describe vacancy formation 
in Ce02, the simplest case to consider is a 2 x 2 x 2 su- 
percell of Ce02 from which a single O atom is removed. 
This supercell has the composition Ce32063 (CeOi.97). 
As this composition does not correspond to any known 
phases of reduced ceria, it is customary to require the 
structural parameters of the relaxed structure to match 
those of Ce02 based on Vergard's rule^, i.e., there is a 
linear relationship between the lattice constant and the 
extent of reduction. As will be discussed in more detail 
in Section |IV C[ most of the DFT work supported the 
standard picture. However, the approximate functionals 
of DFT such as LDA, GGA, and LDA+/7 do not appear 
to describe well the electronic properties of the reduced 
phase d 28 " 30 !. Furthermore, recent calculations using hy- 
brid functionals 31 and DFT+t^MSI do not 

appear to 

support the standard picture. 

As far as we know, a study of the charge distribu- 
tion near bulk O vacancies in crystallographic phases of 
reduced cerium oxide apart from Ce203 (and the cus- 
tomary supercell of DFT calculations, i.e., CeOi.97) has 
not been done. Our study of the charge distribution in 
the well-characterized intermediate phases represents a 
complementary approach to answering the fundamental 
question of where the two electrons left behind when an 
O vacancy forms go. Furthermore, since no in-situ stud- 
ies have been done to identify the precise phases involved 
in ceria-based catalysis, it is also possible that the inter- 
mediate phases may play an important role in the engi- 
neering applications of these materials. 

In this study, we examine the standard picture for the 
oxidation and reduction in ceria described above, consid- 
ering primarily two intermediate phases namely Cen02o 
and Ce70i2- A brief discussion of CeeOn is also in- 
cluded. We discuss this in light of our recent results on 
Ce site valencies obtained by an analysis of observed crys- 
tal structures using the bond valence method 33 . The two 
cases, Cen02o and Ce70i2, were chosen to illustrate the 
evolution of the charge distribution around O vacancies 
that occurs with reduction. Furthermore, apart from the 
extremal structures, namely Ce02 and Ce203, these two 
examples along with Ce^O^ were the only crystal struc- 
tures of the Ce oxides for which we were able to obtain 
highly accurate crystallographic data to enable the type 
of analysis we do here. 

As is discussed later, both Cen02o an d Ce70i2 have 
fluorite-related crystal structures. According to Kang 
and Eyring, Ce70i2 is the limiting intermediate oxide 
for which any further reduction leads to loss of the flu- 
orite structure of Ce02r^. However, the fluorite-related 
CesOs phase was recently observed in the phase diagram 
of the Ce-0 systerrPn In addition to the conventional 
unit cells, we will generally use the Kang and Eyring 
structural principle of modular units to discuss O va- 



cancy ordering in reduced ceria phases ^ 34 * 36 * 37 *. There 
will be oc casional references to Martin's coordination de- 
fect modeP3HHEHl w here this model may provide a better 
conceptual framework. 

This paper is organized as follows: In Section [IT] we 
present the crystal structures of Cen02o and Ce70i2 
both in the conventional way as well as in the represen- 
tation of the Kang-Eyring fluorite-type module theory. 
We show how the Kang-Eyring fluorite-type module the- 
ory helps to understand the structure of O vacancies in 
Cen02o and Ce70i2 in relation to the parent fluorite 
structure. In Section [Hi} we consider the implications of 
our bond valence calculations for the charge distribution 
in the local environment of the O vacancies. Section HVl 
discusses how our results conflict with the standard pic- 
ture of charge localization in reduced ceria phases. We 
then make a prediction of the charge distribution to be 
expected in the CeeOn crystal based on the bond va- 
lence model in Section [Vj The possibility of direct /-/ 
coupling between neighbouring Ce sites in crystals of the 
Ce oxides. We assess the validity of this idea from the 



Harrison method of universal parameters in Section VI 



Conclusions are presented in Section VII 



II. ALTERNATIVE REPRESENTATIONS OF 
THE CRYSTAL STRUCTURES OF CE n 2 o AND 
CE7O12 

A. Coordination defect model 

Martin proposed describing an O vacancy in a Ce oxide 
crystal as a structural entity made up of 4 Ce ions a nd 6 O 
ions around the vacancy site as illustrated in Fig. [|f 4 * 18 l 
The 4 Ce ions form the first coordination shell (a tetra- 
hedron) of the vacancy site while the 6 O atoms form the 
second coordination shell (an octahedron). Compared 
to the perfect structure of Ce02, all 4 Ce ions are dis- 
placed 0.2 A away from the vacant O site while the O 
atoms are shifted 0.3 A towards the vacancy. This struc- 
tural unit has been called the coordination defect and it 
maintains its structural integrity in a crystal of the oxide 
(see Fig. [2J. It is constructed from two types of sub- 
units called octants, one of composition Cei O called the 
S phase and the other incorporating an O vacancy, the 
A phase, of composition Cei. The exploded view in Fig. 
[2] (b) shows these octants and how they are arranged to 
form the coordination defect. This model, in its original 
formulation, does not explicitly account for the changes 
in the valences of the Ce ions and the four Ce ions which 
are nearest neighbours to the O vacancy in the coordina- 
tion defect are implicitly identical. 



B. Modules 

Soon after the elaboration of Martin's coordination 
defect modeP^, Kang and Eyring developed a simpler 



5^. 
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FIG. 2. (a) Coordination defect showing characteristic topol- 
ogy, (b) Exploded view of the coordination defect showing 
the central octant of what is called the "A phase" coordi- 
nated by six MiO octants of what is called the "<5 phase". 
The A phase only differs from the 5 phase in that it has an 
O vacancy as shown. Black circles represent Ce atoms, grey 
(orange in colour version) circles - O atoms and the light grey 
circle inside the A phase is an O vacancy. After Martin™ 



model for constructing the crystal structures of the higher 
rare earth oxides, i.e., oxides from CerOi2 to Ce02 based 
on a different set of structural element^. These so- 
called modules, include the Ce02 unit cell with the rest 
derived from it by creating one or two O vacancies in cer- 
tain prescribed ways resulting in a total of thirteen mod- 
ules. These modules are the fundamental building blocks 
for all the higher oxides under a prescription of rules de- 
fined by the authors 34 . In Fig[3j we show the ten modules 
relevant to Cen02o and Ce70i2 in addition to the unit 
cell of Ce02- Kang and Eyring's method accounted for 
all the higher oxides known at the time and was success- 
ful in predicting the existence of new structures 34 . The 
Kang and Eyring method enables a relatively simple way 
of visualizing vacancy sites which are otherwise difficult 
to decipher when the structure is cast in the conventional 
setting of its space group. This provides considerable fa- 
cility in the analysis of the local environment of the O 
vacancy site. As a result, we will make frequent reference 
to this framework in the rest of this paper. 



u 1 u 2 



u 3 u 4 



D 3 D 4 



FIG. 3. The ten different modules required to build the mod- 
ular unit cells of CenCbo and Ce70i2- The vertices of the 
tetrahedra denote the locations of the Ce ions surrounding an 
O vacancy. Module F is the conventional unit cell of CeC^. 
Except for W± which has two O vacancies along the body 
diagonal of the CeCb unit cell, the rest of the modules have 
only one O vacancy in the positions shown. [Adapted fronP^j 



C. Cen0 2 o 

The conventional unit cell of Cen02o is shown in Fig. 
[4] Note that for all ball and stick figures in this pa- 
per, the sticks between the Ce sites and the O vacancy 
site do not represent chemical bonds. Instead, the sticks 
are only included to show the geometric relationships 
between the atoms. The unit cell of Cen02o consists 
of one formula unit per unit cell with inversion as its 
point group symmetry. It has six Ce and ten O distinct 
sites along with their inversion images with the Ce(l) 
site as the centre of inversion. However, for a descrip- 
tion of the vacancy structure and charge distribution 
in this crystal, it is convenient to view all the lattice 
points in the O sublattice as roughly equivalent and to 
divide the Ce sublattice into two regions. The first re- 
gion consists of the Ce(l) and Ce(2) sites which, in the 
Ce sublattice, occupy the second coordination shells of 
the O vacancies. The remainder of the Ce atoms are the 
first coordination shells of the O vacancies. The crys- 
tal structure of Cen02o can be readily related to the 
fluorite structure of Ce02 from which it is derived by 
reduction. This relationship becomes clearer with the 
Kang- Eyring module representation. Kang and EyrimJ^ 
have shown that the modular unit cell of Cen02o has 
the modular composition 3F, 4D, 4U where the mod- 
ule types are as given in Fig. [3] They also provided 
the modular sequences for the [100] and [010] direc- 
tions in this crystal as FFU 3 U 2 D 3 D 2 FU 4 U 1 D 4 D 1 and 
FU 3 D 3 FU 1 D 1 FU 2 D 2 U 4 D 4 respectively. Using this in- 
formation along with their module juxtaposition rules, we 
derived the modular sequence for the [001] direction to 
be D^FFFU^U^D+Ds. This then enabled us to 
construct the 3D modular unit cell of Cen02o shown in 
Fig. [5j We note that this unit cell is not unique. There is 
only one type of O vacancy in the crystal of Cen02o and 
if we only consider the Ce sublattice, each O vacancy has 
the four Ce(3)-Ce(6) sites forming its first coordination 
shell. These four Ce sites form a distorted tetrahedron. 



Ce{l),Ce{2) 
Ce{3)-Ce{6) 
O 

O vacancy 



FIG. 4. The unit cell of CenCbo whose space group is PI. 
There are six distinct Ce sites in this crystal. However, to fa- 
cilitate the discussion of vacancy structures and charge distri- 
bution, the Ce sublattice has been partitioned into two types 
of lattice points as shown in the legend. The first set of lat- 
tice points consists of the Ce(l) and Ce(2) sites which, if we 
only consider the Ce sublattice, occupy the second coordi- 
nation shells of the O vacancies. The second set consists of 
the remainder of the Ce sites and these form the first coor- 
dination shell of the O vacancies in the Ce sublattice of the 
crystal. There are ten distinct O sites in the crystal but for 
our purposes, it is convenient to treat all the O sites in the O 
sublattice as roughly equivalent as indicated in the legend. 




(a) 



(b) 



(c) 



FIG. 5. The fluorite-type modular unit cell of Cen02o de- 
rived using the Kang and Eyring rulesP. In (a), the modular 
unit cell is constructed from transparent modules to show the 
location of the O vacancy site as indicated by the shaded 
tetrahedra with the module types as identified in (c) and 
shown in detail in Fig. [3] A simplified block representa- 
tion of the modular unit cell given in (a) is shown in (b). The 
exploded view with all the module types identified is given 
in (c). The modular composition is: top layer - U 1 , F, U 4 ; 
middle layer - Di, D±, F, F, U 3 ; bottom layer - D 2 , D 3 , U 2 . 
These figures make clear the distribution of the O vacancies 
within the crystal in relation to the original fluorite structure. 
There are large contiguous regions, modules F, where there 
are no O vacancies. 



FIG. 6. The proposed vacancy cluster of CenCbo consisting 
of two O vacancies which are nearest neighbours and related 
to each other by inversion symmetry centred on the Ce(l) 
site. The various Ce sites are labelled as shown. The white 
and light grey (teal) circles represent O vacancies and various 
O sites respectively. 



The O vac - Ce(i) distances are shown in Table [T] where 
i = 3,4,5,6. The second coordination shell consists of 
twelve Ce sites namely, two Ce(l), four Ce(2), two Ce(3), 
two Ce(4) and one of each of the Ce(5) and Ce(6) sites 
as shown in Table U Of these twelve Ce sites in the O va- 
cancy's second coordination shell, only half a Ce(l) site 
and one Ce(2) site can be assigned to the particular O 
vacancy site. Thus, in total, an O vacancy in Cen02o is 
proportionately associated with five and half Ce sites. It 
is convenient to work with an integral number of Ce sites 
when studying the charge distribution in the local envi- 
ronment of an O vacancy. In order to do this, we propose 
to partition the lattice so that we consider a vacancy clus- 
ter consisting of two neighbouring O vacancies which can 
then be associated with eleven Ce sites. Since the choice 
of the Ce(l) site is the most difficult to establish, we pick 
a Ce(l) site and construct the vacancy cluster around it. 
A Ce(l) site has four closest O vacancies all in the third 
coordination shell. Two of the O vacancies are located 
4.558A away from the Ce(l) lattice point with the other 
two at 4.571A. Since the Ce(l) site is a centre of in- 
version, each of these pairs of O vacancies are related by 
inversion symmetry. We choose the pair with the shortest 
distance from the Ce(l) site to define the vacancy clus- 
ter consisting of two O vacancies. The cluster obtained 
by this procedure is illustrated in Fig. [6j This cluster 
should not be confused with the divacancy of CerOi2 to 
be discussed later. The divacancy of Ce70i2 is a unique 
structural entity of the crystal arising from long-range or- 
dering of the vacancies whereas the cluster we define here 
is simply an analytic convenience. It cannot be strictly 
considered to be specified in a unique way since the pair 
of O vacancies we chose is only O.OlA closer to the Ce(l) 
compared to the other. On the O vacancy sublattice, 
each O vacancy has four nearest neighbours 6.134A away 
for the first pair and 6. 138 A for the second. There is 
no difficulty in identifying the Ce(2) site associated with 
each O vacancy since there is only one closest (4.495A) 
Ce(2) site. The distances are calculated from the crys- 
tallographic data of Kummerle and Hegei 40 . Zhang et 
al^ suggested another useful way of viewing the crys- 
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TABLE I. The first- and second coordination shells of the O 
vacancy in the Ce sublattice of CenCbo and some selected dis- 
tances between the sites. These Ce-Ce distances show that the 
average intra-tetrahedral separation between Ce sites (4. 2 A) 
is larger than the inter-tetrahedral separation (3.9A). 

Ce site Distance from Direct Ce-Ce Direct Ce-Ce 

O vacancy, A type distance, A 





First Coordination 




Ce(6) 


2.5106 


Ce(6)-Ce(3) 


4.169 


Ce(3) 


2.5674 


Ce(6)-Ce(4) 


4.194 


Ce(4) 


2.5817 


Ce(6)-Ce(5) 


4.207 


Ce(5) 


2.6185 


Ce(3)-Ce(4) 


4.219 






Ce(3)-Ce(5) 


4.144 






Ce(4)-Ce(5) 


4.234 




Second 


Coordination 




Ce(5) 


4.3184 


Ce(l)-Ce(5) 


3.810 


Ce(4) 


4.4140 


Ce(l)-Ce(4) 


6.904 


Ce(3) 


4.4164 


Ce(l)-Ce(3) 


3.852 


Ce(3) 


4.4341 






Ce(4) 


4.4390 






Ce(6) 


4.4438 


Ce(l)-Ce(6) 


5.561 


Ce(2) 


4.4952 


Ce(l)-Ce(2) 


7.715 


Ce(2) 


4.5265 


Ce(2)-Ce(3) 


3.873 


Ce(2) 


4.5340 


Ce(2)-Ce(4) 


3.861 


Ce(l) 


4.5579 


Ce(2)-Ce(5) 


6.860 


Ce(l) 


4.5705 


Ce(2)-Ce(6) 


5.573 


Ce(2) 


4.6007 







tal structure of Tbn02o which is also directly applicable 
to Cen02o because these compounds are isostructural. 
Their perspective is based on the observation that the O 
vacancies are more or less uniformly distributed within 
the crystal structure with no evidence of vacancy pair- 
ing. The vacant O sites are distributed in a way such 
that the separation between the defects is maximized and 
thus, the repulsive interactions are reduced. The Ce ions 
in the first coordination shell relax outward in a more 
or less isotropic way (Ce-F ' ranges 2.51 — 2.62A, see 
Table l|). As the vacant Ce tetrahedra expand, the re- 
gions between them experience substantial compression 
as indicated by the difference in the distances between 
the intra-tetrahedral ions (4. 2 A) and those of the inter- 
tetrahedral ions (3. 9 A). These distances are given in Ta- 
ble [TJ In this table, the intra-tetrahedral ion distances 
are given as the direct Ce-Ce distances in the first co- 
ordination shell while the inter-tetrahedral ion distances 
refer to these distances in the second coordination shell: 
Ce(l)-Ce(5), Ce(l)-Ce(3); Ce(2)-Ce(3) and Ce(2)-Ce(4). 
The inter-tetrahedral Ce ions have a coordination num- 
ber of 8 as in the parent fiuorite structure. In discussing 
the charge distribution in the local environment of the 
vacancy, we will notice that the charge is distributed in 
a manner which is distinctly different between the intra- 
tetrahedral and inter-tetrahedral Ce ions. 




FIG. 7. The unit cell of CerOi2 showing the two distinct Ce 
sites: Ce(l) - black and Ce(2) - dark grey. The 0(1) sites are 
shown in grey (teal) with the 0(2) light grey (orange). The 
Ce(l) and Ce(2) sites are of 5*6 and i symmetries, respectively. 
The O vacancy sites are shown in white. A characteristic 
structural unit of the CerOi2 crystal is the divacancy, which 
consists of two O vacancy sites connected by an $6 Ce site 
between them. 



W, 3 at centre 




(a) 



(b) 



FIG. 8. (a) The modular unit cell of CerOi2 showing the de- 
tails of the seven modules used to construct it. Each tetrahe- 
dron represents the four Ce atoms around the oxygen vacancy, 
(b) A simplified model of the modular unit cell. Adapted from 
Kang and Eyring 34 . 



D. Ce 7 Oi2 

Fig. [7] shows the unit cell of Ce70i2 in the conven- 
tional rhombohedral setting. In the Kang-Eyring frame- 
work, the modular composition of this crystal is D 2l 
L> 4 , W?, U 1 , U 2 and U^. The modular unit cell con- 
structed from this basis is depicted in Fig. [8j It is 
readily apparent from Fig. [8] where the O vacancy sites 
are located. There are two types of Ce sites, Ce(l) and 
Ce(2), in the crystal lattice of Ce rOi2 which are of Sq 
and i symmetry respectivel y 40 * 42 *. These distinct sites 
are indicated in the conventional unit cell of Fig. [7[ Al- 
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though Kang and Eyring did not explicitly identify the 
distinct sites in their modular unit cell according to their 
symmetries, it is apparent that by translating their unit 
cell in three dimensions, two distinct Ce sites are ob- 
tained; the $6 site with two nearest neighbour O vacan- 
cies around it and the i site with only one O vacancy. 
The 5e site has been called the divacancy site^ and it 
forms a shared corner between two coordination tetrahe- 
dra of the O vacancies as illustrated in Fig. [TT] We list 
here some interatomic distances that we will use in our 
analysis of the divacancy cluster: Ce(2)-Ce(2) = 4.09A, 
Ce(l)-Ce(2) = 4.11A, Ce(l)-0(2) = 2.19A and for Ce(2)- 
0(1), there are two relevant bond lengths, namely 2. 23 A 
and 2.32A. The Ce-Ce distances are direct distances be- 
tween these sites, not distances along the Ce-0 bonds 
connecting them. It appears that the Ce sublattice has 
not changed much relative to Ce02- 

It is helpful to visualize vacancy ordering in Ce70i2 
by viewing the divacancy as a structural unit of the O 
vacancies in this crystal. This clarifies why the ratio of 
Sq to i sites is 1 : 6; all the Ce sites which are nearest 
neighbours to the two O vacancies forming a divacancy 
are of Sq symmetry. As there are seven Ce atoms in 
the formula unit, the rest of the Ce atoms (six) must 
be the Ce(2) sites of i symmetry. The formation of the 
divacancy structures as Ce02 is reduced to Ce70i2 is 
associated with some relaxation of the crystal lattice al- 
though the Ce sublattice remains somewhat invariant!^. 
We now examine this relaxation process which is illus- 
trated in Fig. [9] and based on Kummerle and Heger's 
neutron diffraction datsP^. The bond lengths in Ce02 are 
2.434A and 2.706A for the Ce-0 and 0-0 bonds respec- 
tively. When a divacancy is formed, the bonds within 
the divacancy relax and diffraction data show that in 
Ce 7 0i2 the 'bond' lengths are 40 : 2.505A for 0(1)-V£, 
2.364A for 0(2)-^, 2.418A for Ce(l)-V ' and 2.543A for 
Ce(2)-V G *. This structural data indicates that compared 
to Ce02, the bonds within the divacancy have changed 
as follows: the 0(1)-Vq and 0(2)-Vq 'bonds' contract by 
0.20A and 0.34A respectively while the Ce(2)-V ' 'bond' 
gets longer by 0.109A. The Ce(l)-Vo 'bond' is prac- 
tically unchanged, showing only a small contraction of 
about 0.02A. Thus, we see that, just as was the case for 
Cen02o, the O atoms are attracted to the vacancy site 
while the Ce(2) atoms are repelled away so that the O 
atoms now become the nearest neighbours of the vacancy 
site. This appears to be a general feature of vacancies 
in rare earth higher oxides^^U. We note that the six 
O atoms around the vacancy relax by different amounts 
with the 0(2) atoms being more strongly attracted to the 
vacancy site. We also note that according to the coordi- 
nation defect model, the crystal structure of CerOi2 has 
the highest packing density of the coordination defect. 
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FIG. 9. The relaxation of atoms forming part of the divacancy 
in CerOi2. The different lattice sites are indicated by the 
colours as follows: Ce(l) - black, Ce(2) - dark grey, 0(1) 
- grey (teal), 0(2) - light grey (orange) and the O vacancy 
(Vo) - white. The arrows indicate the direction in which 
the respective atoms relax relative to the O vacancy site and 
the amounts by which the 'bonds' relax are shown in the 
legend. This figure was produced from the crystal structure 
determined by neutron scattering 40 . 



III. EXCESS ELECTRONS DELOCALISE AWAY 
FROM THE OXYGEN VACANCY 

Having discussed the O vacancy ordering in Cen02o 
and CerOi2 as well as proposing some vacancy clusters 
suitable for the analysis of charge distribution in these 
crystals in the previous section, we now present the re- 
sults of the bond valence calculations, I, in these clusters. 



Cen0 2 



Fig. 



10 shows the bond valence sums on the various 
Ce sites in the vacancy cluster of Cen02o which was in- 
troduced earlier, see Fig. [6] Clearly the charge from the 
vacancy does not localize on the Ce sites closest to the 
vacancy but rather, delocalizes onto the Ce(l) and Ce(2) 
sites. This contradicts the standard picture. The results 
of Fig. 10 can be compared with those of Tbn02o with 



which it is isostructural. Based on a simple consideration 
of the crystal radii of the cations^, Zhang et at™ sug- 
gested the following valencies for the cations: Tb(l) and 
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FIG. 10. Charges of the different Ce ions in the vacancy 
cluster of CenCbo obtained from the bond valence model. 
The vacancy cluster consists of two vacancies connected by 
the Ce(l) site which is also a centre of inversion for the cluster. 
The four excess electrons that arise due to the removal of the 
two oxygen atoms do not localise on the cerium atoms nearest 
the vacancies. This is contrary to the standard picture. 



Tb(2) - +3, Tb(3) and Tb(4) - +3.75, Tb(5) and Tb(6) 
- +4. Here the numbering of the sites is the same as 
for Cen02o given in Fig. [6] Their assigment of the site 
valences qualitatively agrees well with our bond valence 
results for Ce n 2 o. The Tb 3+ and Tb 4+ ions are / 8 and 
f configurations respectively and can be compared to 
the corresponding f 1 and f° configurations respectively 
for the Ce ions. They suggested the possibility of fast 
electron transfer between the neighbouring Tb(2) sites 
so that, instead of two Tb(2) 3+ and one Tb(2) 4+ instan- 
taneous states, an average is obtained from the neutron 
diffraction data. It was suggested by the authors that the 
thermal ellipsoids of the atoms in Ce70i2 which are much 
higher than in Ce02 could be evidence of the dynamic 
disorder arising from the fast electron transfer. 



B. Ce 7 Oi2 

The charge distribution in the divacancy of Ce70i2 is 
shown in Fig. [TT] This figure clearly shows that, apart 
from the Ce(l) site which is of a different local site sym- 
metry compared to the Ce(2) sites, the charge delocal- 
izes over all three of the Ce atoms closest to the vacancy. 
Again, this contradicts the standard picture, which has 
the electrons completely localise on two of the neighbour- 
ing Ce atoms. 



IV. FAILURE OF THE STANDARD PICTURE 
OF CHARGE LOCALISATION NEAR A 
VACANCY 

We now discuss the above results comparing them to 
the standard picture and propose an alternative view of 
how the charge redistributes itself following O vacancy 
formation in ceria as deduced from the bond valence 
model. As already mentioned, the conventional descrip- 
tion of the electronic processes involved during O vacancy 
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FIG. 11. The divacancy of Ce70i2 showing the valences of 
the Ce atoms calculated from the bond valence model, I. The 
charge delocalizes on the three Ce(l) sites of triclinic symme- 
try which now have a valence of +3.21. 



formation requires that the excess charge at the vacant O 
site be shared between two of the four Ce ions which form 
the first coordination shell of the vacancy. This means 
that the tetrahedron formed by the four Ce ions around 
the O vacancy consists of two Ce 3+ and two Ce 4+ ions. 
This description is based on the ionic model for these 
oxides although Kang and Eyring do indicate that there 
may be some covalent character in these oxides 34 . 



A. Cen0 2 o 

If we assume that the four electrons which are left 
when the two O vacancies are formed are localized within 
the cluster, then we can describe how the charge redis- 
tributes itself. In fact, from the bond valence results, 
which also showed mixed valence in Ce02, the valence 
of Ce in Ce02 is +3.73 and the corresponding valence of 
O is —1.87. A description of vacancy formation in Ce02 
to form Cen02o can then be illustrated in the bond va- 
lence model as shown in Fig. 



12 Thus, according to 



the mixed valence description of the Ce02 crystal from 
the bond valence model (see Fig. 12 'a)), when an O 
vacancy is created in Ce02, a charge of —1.87 is left be- 
hind as opposed to a charge of —2.0 in the ionic picture. 
This then means that, if we consider the vacancy clus- 
ter of Fig. 10 and equivalently, Fig. 12 'b), the creation 
of the two O vacancies leaves a total charge of —3.74 in 
the cluster. From Fig. [lO} we can see that compared 



to their original valences in Ce02, the valences of Ce(5) 
and Ce(6) are virtually unchanged. Thus, considering 
the change in the valences of the remainder of the Ce 
sites, the total charge accummulation on the Ce sites in 
the cluster is —2.87. There is therefore a discrepancy 
of ~ — 1 between the total extra charge in the Cen02o 
crystal inferred from the bonding in Ce02 (i.e. —3.74) 
and that directly calculated for this crystal by the bond 
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+3.73 



+3.58 



+3,06 



FIG. 12. A schematic to illustrate the formation of Cen02o 
from CeC>2 through vacancy formation in the O sublattice. In 
(a), a part of the CeCb lattice, which, after removing the two 
O atoms (light grey (orange) circles), gives a cluster corre- 
sponding to the one given in Fig. [6] Also note that compared 
to Fig. [6] here we have made the further simplification of re- 
placing all the Ce-O-Ce bonds with direct Ce-Ce bonds except 
for the two O atoms which are subsequently removed to create 
vacancies. For these two atoms, no bonding detail is given in 
the schematic except that we place each of them at the centre 
of a regular Ce tetrahedron reflecting their coordination en- 
vironment in Ce02- The numbers shown at the atomic sites 
are the valences of the respective atoms as determined by the 
bond valence methoo™, and Fig. 
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In (b), the cluster of (a) 
after reduction to create the two O vacancies (empty circles) . 
The colours in this figure denote the following: black circles - 
Ce sites which are not nearest neighbours to any O vacancies, 
grey circles - Ce atoms which are the nearest neighbours of 
at least one O vacancy, light grey (orange) circles - O atoms 
and empty circles - O vacancies. 



valence calculations (—2.87). This discrepancy is entirely 
consistent with the accuracy of the bond valence model 
reported in I although we note that the accounting of 
charge we have done here from Ce02 to Cen02o is very 
simplistic. It ignores the many-body effects which may 
follow the injection of extra charge at the Ce sites when 
the O vacancies are created. Such an extra charge may 
change the degree of mixing between the Ce 4/- and the 
O 2p states in manner that may not be related in a sim- 
ple way to their state in Ce02- It may happen that, for 
instance, due to a change in the hybridization of Ce 4/- 
and the O 2p states, the extra charge may not be confined 
to the Ce 4/ as we assumed in calculating the discrep- 
ancy in the total charge, but rather, some of the charge 




FIG. 13. Schematic illustrating the charge distribution fol- 
lowing divacancy formation, based on a bond valence sum 
analysis of the crystal structures of CeCb and Ce70i2- (a) 
Part of the CeCb lattice showing two neighbouring O atoms 
(light grey (orange) circles), each inside a tetrahedron of Ce 
atoms (black circles). The numbers shown are the bond va- 
lence sums for each of the atoms. When the two neigbouring 
O atoms are removed due to the reduction of CeCb, the diva- 
cancy of CerOi2 is formed as shown in (b). The bond valence 
sums on the Ce sites of triclinic symmetry (grey circles) have 
decreased from +3.73 in Ce02 to +3.21 reflecting the presence 
of the additional charge left behind when the two O atoms are 
removed. The bond valence sum on the 5*6 Ce site is virtually 
unchanged. 



may go into the O valence band. It is expected that the 
results in Fig. [lOjwould capture these many-body effects. 
An important point from the results of Fig. [To] is that 
-1.99 (~ 70%) of the total charge of -2.87 goes to the 
three Ce atoms in the second coordination shell with the 
remainder being shared among the eight Ce atoms in the 
first coordination shell. 

If one just considers the cation sublattice, it is evi- 
dent that the charge accumulates in regions of high strain 
(inter-tetrahedral Ce ions) where the Ce-Ce distances are 
slightly shorter compared to the intra-tetrahedral case. 
This led some authors to consider the possibility that 
there may be some direct /-/ hopping in these dense 
regions 41 . As we will discuss in Section VI for the case 
of CeyOi2, direct /-/ hopping appears unlikely in these 
oxides. 



B. Ce 7 Oi2 

Based on the results of Fig. [TT] we can now discuss the 
electronic features of the local environment of the O va- 
cancy in Ce70i2- The formation of the Ce70i2 divacancy 
from the Ce02 crystal lattice according to the BVM is 
schematically illustrated in Fig. 13 'a) and (b). It can 



be seen that two corner-sharing Ce tetrahedra of Ce02 
lose an O atom each from the centre of the tetrahedron. 
Our bond valence calculations show that the valences of 
the Ce(l) and Ce(2) sites in CerOi2 are +3.67 and +3.21 
respectively. The valence at the Ce(l) site is comparable 
to the value of +3.73 calculated for Ce in Ce02 using a 
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bond valence method analysis of the crystal structure. If 
we use the Kang and Eyring modular unit cell with its 8 
O vacancies, 4 Ce(l) and 24 Ce(2) sites and further as- 
sume that we start from mixed-valence Ce02 where the 
valences are +3.73 and —1.87 for the Ce and O sites 
respectively, we can get some estimate of how the excess 
charge on the O vacancy sites is distributed in the mod- 
ular unit cell. The valence of the Ce(l) sites does not 
change much from its value before the O vacancy is cre- 
ated. Thus, virtually all the charge left at the vacancy 
site by an O atom is evenly distributed within the first 
coordination shell of the vacancy site among the (three) 
Ce(2) sites. On this basis and only considering a sim- 
ple analysis using the bond valence results for Ce02, the 
valence at the Ce(2) sites is then readily calculated to 
be +3.11 (i.e., 3.73 — 1.87/3) which is comparable to the 
independently calculated value of +3.21 within the error 
bounds (±0.1) of the method. That the excess charge is 
not distributed evenly among all four of the Ce atoms in 
the first coordination shell contradicts the standard pic- 
ture. If the charge was evenly distributed among all four 
Ce atoms, the valences would be +2.79 and +3.26 for the 
Sq and i sites respectively. 

The above picture from the BVM can be contrasted 
with standard picture description of the same process as 
illustrated in Fig. 13 'c) and (d). Here, there are two 



ways one could assign the valences of the Ce sites in the 
divacancy. One way is to assign valences of +4 to three 
of the i sites, two on the first tetrahedron and one on the 
second. The remainder of the Ce sites are then assigned 
the valence of +3 as illustrated in Fig. [l3^d)(i). An 
alternative arrangement of the charges is illustrated in 
Fig. 13^d)(ii). Here, the three +4 valences are assigned 
one each to one of the i sites from the first and second 
tetrahedra with the third being assigned to the Sq site. 



C. Comparison with calculations based on Density 
Functional Theory 

Most DFT calculations (performed in the supercell 
of composition Ce32063) characterizing the location and 
character of the excess charge near bulk O vacancies con- 
cluded that the charge localize s on two of the nearest 
neighbour Ce ions 12 ' 15 ' 16 ' 26 ^^. This is what we have 
called the standard picture. It appears that this picture 
was first brought into question by Castleton et alW^ in 
a paper where the question of charge localization was 
extensively explored. They concluded that it wasn't pos- 
sible to fully localize the Ce 4/ electrons on these Ce sites 
in the DFT+£7 framework while at the same time pre- 
serving a correct description for all the other electrons. 
We are not aware of any subsequent DFT+Z7 reports di- 
rectly addressing this question in the conventional super- 
cell. The more recent work used a periodic electrostatic 
embedded cluster method (PEECM)pA This work found 
that the charge preferred to localize well away from the 
O vacancy, being on Ce sites of the third coordination 



shell. As this shell coincided with the boundary of the 
quantum mechanical part of the cluster used in the cal- 
culation, the authors raised the possibility that the result 
may be an artifact of the method. 

We also note that Burow et a/P^ obtained three differ- 
ent Ce-Ce distances of 4.07, 4.12 and 4. 20 A in the first co- 
ordination shell. These distances are shorter than those 
reported in Table [I] for the vacancy cluster of Cen02o- 
However, since the structures in question are different, 
one should be cautious about such comparisons. 

Our results are also in agreement with some recent re- 
ports on the location and localization of th e exce ss elec- 
trons for surface and subsurface O vacancie d 31 * 32 1 . These 
authors argued that the localization of charge on next 
nearest neighbour Ce sites was mainly controlled by the 
lattice relaxation due to electrostatics as reflected in the 
lowering of the Madelung potential. However, it should 
be remarked that a bulk O vacancy is in a different co- 
ordination environment relative to either a surface or 
subsurface O vacancy. It is therefore possible that the 
comparison between our results for bulk O vacancies 
and those reported for surface/subsurface O vacancies 
may not be simple. Indeed, some authors have reported 
differences in the location of the charge between sur- 
face/subsurface and bulk O vacancies^. Nevertheless, 
it is interesting to note that, until these more recent re- 
sults, the standard picture was considere d applic able to 
surface and subsurface O vacancies as well 1 ^^^. 



V. PREDICTED CHARGE DISTRIBUTION IN 

CEeOn 

We make a few remarks about a possible charge dis- 
tribution around the O vacancies in this crystal whose 
unit cell is shown in Fig. 



14 Complete crystal struc- 



ture data is not available but we make predictions of the 
charge distribution to be expected in this crystal from 
a generalization of the bond valence results of Ce70i2 
and Cen02o- The proposed vacancy cluster for study- 
ing the charge distribution in the local environment of 



the oxygen vacancy in this crystal is shown in Fig. 15 
We did not find a full crystallographic characterization 
of this crystal and therefore it is not possible to perform 
accurate bond valence calculations for it. The results re- 
ported in I for this crystal were based on the positional 
parameters of Pr 6 On. If one were to accept these re- 
sults (Table I of I), then since Ce(l), Ce(2), Ce(3) and 
Ce(4) coordinate the O vacancy with Ce(5) and Ce(6) in 
the second coordination shell for the Ce sublattice, we 
see that the charge distribution does not quite follow the 
pattern one would expect from the results obtained for 
Cen02o and Ce70i2- The overall charge in the vacancy 
cluster is larger by about ~ — 1.2 compared to what one 
would expect starting from mixed valence Ce02- This 
discrepancy is twice what we found in the case of Cen02o 
if the discrepancies are normalized to the total number 
of Ce sites in each vacancy cluster. As a result, the three 
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Ce(3) 



£Ce(l)-Ce(4) ^Ce(5),Ce{6) O 



O vacancy 



FIG. 14. The unit cell of CeeOn which consists of four for- 
mula units. The Bravais lattice is monoclinic and of space 
group P2\/c. 




Ce{5) 

Bond Valence Estimates 
Ce(l):+3.29 Ce(4):+3,70 
Ce(2):+3.47 Ce(5):+2.93 
Ce(3):+3.09 Ce(6):+3.14 



Predicted Valences 
Ce(l)-Ce(4):+3.7 
Ce(5),Ce(6): +3.0 



Ce sites Ce(l), Ce(2) and Ce(3) which one would expect 
to have valences close to +3.7 have much lower valences. 
We do not believe that this is the correct charge distri- 
bution in the local environment of the O vacancy for this 
crystal. Instead, we prefer to extrapolate the argument 
which emerges from the bond valence results of Cen02o 
and CerOi2. Doing so leads to the prediction that the 
correct charge distribution for this crystal should be as 
illustrated in Fig. [15] where all the charge goes to the 
two Ce sites in the second coordination shell i.e. Ce(5) 
and Ce(6). Of the known crystal structures of the re- 
duced higher oxides of ceria, CeeOn is the closest to the 
customary supercell used in DFT in that it is the least 
reduced. The prediction we have made here would sug- 
gest that, except in an artificial way, it should not be 
possible to localize the electrons from the O vacancy site 
in the first coordination shell for the customary supercell 
used in DFT for studying this problem. An experiment 
giving the full crystallographic data for this phase would 
help resolve this question. 



VI. DELOCALISATION VIA f-p 
HYBRIDISATION 



The possibility of the direct hopping of electrons be- 
tween /-orbitals on neighbouring Ce sites has been sug- 
gested to explain the charge distribution in reduced ceria 
phases 41 . We briefly explore this suggestion considering 
the example of Ce70i2- There are two immediate ques- 
tions one may like to answer about the charge distribu- 
tion obtained for this phase from the bond valence calcu- 
lations. First, how does an /-electron delocalize among 
the Ce(2) sites? Is this through direct /-/ hopping or 
through the more indirect two-step process of f-p hop- 
ping which involves an O site bonded to both Ce sites? 
Second, why does this derealization of the charge not 



FIG. 15. The proposed vacancy cluster of CeeOn consisting 
of the six Ce sites as shown. The sites Ce(l) - Ce(4) are in the 
first coordination shell of the O vacancy and are coordinated 
to 7 O atoms. The Ce(5) and Ce(6) sites are in the second 
coordination shell and have the full 8-coordination of the Ce 
site in Ce02- Bond valence sums are listed for each of these Ce 
sites as estimated from the positional parameters of Pr 6 On 
and taken from I. We also show the site valences we predict 
from a generalization of the bond valence results of CerOi2 
and Cen02o- 



extend to include the Ce(l) which is also in the same 
first coordination shell of the O vacancy? We are not in 
a position to give a detailed answer to these questions in 
the present paper but only some preliminary indications 
based on Harrison's method of universal parameters^. 
We do not believe that the derealization of the /-level 
charge between the Ce(2) sites is a result of direct /-/ 
coupling between these sites. 

We now consider the two relevant cases, i.e., direct /-/ 
hopping between neighbouring Ce sites and the indirect 
f-p hopping which involves an O site between the Ce 
sites in question. In each of these cases, we consider 
two situations: electron hopping between two Ce(2) sites 
which delocalizes the charge between these sites and the 
electron hopping between a Ce(2) and a Ce(l) site which 
could also delocalize the charge between these sites. 

For the case of direct /-/ coupling, the relevant dis- 
tances are: Ce(l)-Ce(2) - 4.11A. These distances are 
quite comparable and they give the sa me dir ect /-/ hop- 
ping matrix element, tff, be estimate d 49 * 50 * to be about 
0.01 eV. This matrix element is the most favourable of 
all the four nonvanishing /-/ matrix elements and rep- 
resents a a-<j interaction between the two orbit als. This 
result would suggest that, if direct /-/ coupling were 
the dominant mechanism of charge derealization, then 
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charge would delocalize over all four Ce atoms which are 
nearest neighbours to the O vacancy. 

In the case of f-p hopping, we again, consider two 
cases: hopping along the Ce(l)-0(2)-Ce(2) bonds which 
we designate '0(2) hopping' and the 'O(l) hopping which 
occurs along Ce(2)-0(1)-Ce(2). For 0(2) hopping, the 
respective matrix elements, tf p , for electron hoppin g be- 
tween Ce(l) and 0(2) as well as 0(2) and Ce(2) ar<PEU 
0.67 and 0.46 eV. These matrix elements refer to a a- 
<j interaction between the two orbitals which are the 
more favourable of the two nonvanishing matrix ele- 
ments. Since, for f-p hopping to couple / states between 
two Ce sites, a two-step process involving the hopping of 
an /-electron from one Ce(2) site to an O 2p level fol- 
lowed by the hopping of an electron from the O 2p to the 
4/-level of the other Ce(2) site is required, the overall 
matrix element for 0(2) hopping, £ e //> becomes 0.15 eV 
calculated from Eq. 



Where we have assumed that the energy gap between the 
Ce 4/- and the O 2p levels, e f -e p , is - 2 eV 30 . The cor- 
responding matrix elements for 0(1) hopping 0.50, 0.61 
and 0.15 eV respectively. We see that the overall matrix 
element, t e ff, is the same in the two cases. Again, this 
would suggest that the charge should completely delocal- 
ize over all the Ce sites in the divacancy which contradicts 
the results of the bond valence method. Thus, we con- 
clude that the Harrison parameters cannot fully explain 
the charge derealization in the Ce70i2 divacancy. How- 
ever, they appear to exclude the possibility of direct /-/ 
hopping since t e ff is an order of magnitude larger than 
tff. These approximate results give some indication of 
why the charge derealization may not extend to include 
the Ce(l) site. An alternative way of looking at the in- 
direct /-/ coupling is to view it as a hybridization of 
the 4/ states at the two Ce(2) sites and the O 2p states 
at the O atom connecting them through Ce-0 bonds. 
This indirect coupling is analogous to the superexchange 
mechanism first discussed by McConnelP^. 

VII. CONCLUSIONS AND FUTURE 
DIRECTIONS 

The goal of this review was to answer the question: 
when an oxygen atom is removed from bulk cerium ox- 
ide where do the two excess electrons go? The approach 
taken was to consider high resolution crystal structures of 
Ce n 02o and Ce70i2 and see how they could be viewed 
as ordered arrays of oxygen vacancies in an underlying 
Ce02 crystal. The charge distribution in the local envi- 
ronments of the O vacancies can then be deduced from 
the bond valence model. An important conclusion is that 
the results are incompatible with the widely accepted 
standard picture of charge localization on two cerium 
ions next to the vacancy. Instead, we found that the 



charge distributes itself predominantly in the second co- 
ordination shell of cerium ions. Furthermore, one excess 
electron can be delocalised over more than one cerium 
ion. 

Our conclusions concerning the charge distribution 
near oxygen vacancies are significant for several reasons. 
First, they contradict many (but not all) atomistic sim- 
ulations based on density functional theory. Second, the 
actual charge distribution around the defect has impor- 
tant implications for the other questions we posed at the 
beginning of this review. The charge distribution has a 
significant effect on the relative stability of surface and 
subsurface vacancies^. Also, if the charge around oxygen 
surface and subsurface vacancies is not simply localised 
on Ce ions next to the vacancy this could change our 
understanding of the catalytic activity of these surfaces 
since it has been claimed or assumed that it is associ- 
ated with Ce 3+ ions at the surface.^ 5 - 3 -! The charge dis- 
tribution around the vacancies has implications for elec- 
tronic and ionic conduction, a subject we have discussed 
elsewhere 54 . 

Since our determination of the charge distribution 
around vacancies is indirect it will be important to see 
whether the same conclusions are obtained on experi- 
ments on these higher order oxides with complementary 
probes such as Electron Spin Resonance (ESR), nuclear 
magnetic resonance, and photoemission spectroscopies. 
For example, ESR spectra should quite be distinct for 
the standard picture of electrons localised onto Ce 3 + ions 
and our alternative picture of partial delocalisation over 
the second co-ordination shell and unique charge (and 
spin) distributions associated which single vacancies and 
di vacancies. 

A variety of spectroscopic probes have shown that the 
formation of oxygen vacancies is also associated with next 
electronic states located in energy between the valence 
band of nominally "oxygen 2p" states and the "localised 
4f states ! 55 ^ 56 ! DFT-based calculations which do not suf- 
ficiently include the effect of electronic correlations fail 
to produce these states PS Generally the energy of defect 
states is correlated with the extent of electron localisation 
around defects. Hence, we suggest that spectroscopic 
studies of ordered phases of higher order oxides could be 
a fruitful approach to characterising the electronic prop- 
erties of defects and complementary to the structural ap- 
proach we have reviewed. We expect that the energies of 
the electronic states associated with each of the vacancies 
we have considered will be different and correlated with 
the extent of charge delocalisation that we find. Such 
experimental results will provide significant constraints 
on theories. 

This review only considered bulk cerium oxides. How- 
ever, we believe our approach to characterising electronic 
properties of oxygen vacancies in terms of a bond valence 
sum analysis and the physical insights gained should also 
be fruitful in the study of surface defects and to other 
widely studied oxidessuch as those of titanium, hafnium, 
vanadium, and zirconium. 
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